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CN . I. INTRODUCTION 

o 
o 

Double beta decay is of a major importance both for particle and nuclear physics. In particular, the lepton number 
violating neutrinoless mode of double beta decay (Oz^/3/3-decay) , forbidden in the Standard Model(SM), touches upon 
: one of the key questions of particle physics: if massive neutrinos are Majorana or Dirac particles. As is known this 
process occurs if and only if neutrinos are Majorana particles. Therefore any experimental observation of Oi^/3/?-decay 
would be an unambiguous indication of the Majorana nature of neutrinos even if neutrino contribution to its rate is 
' negligible in comparison with some other mechanisms related to physics beyond the SM Jj . It is commonly believed 
'-H \ that the future double beta decay experiments have realistic chances to observe 0i//3/9-decay and, thus, to establish 
^ 1 the Majorana nature of neutrinos, provide an important information on neutrino mass matrix and on the lepton 
(— ( ' number violating parameters of physics beyond the SM. The accuracy of the determination of these quantities from 
• • I the measurements of the Oi^/3/3-decay half-life will be essentially related to our knowledge of the corresponding nuclear 
_ ■ matrix elements. Therefore, the development of reliable models for calculation of these matrix elements becomes 
, more and more challenging problem of nuclear physics. One of the guiding principles for construction^f such models 
^ ' consists in testing their predictions for the SM allowed two neutrino double beta decay (2t'^-decay) '2', which has 
9^. \ been directly observed in ten nuclides ^J, including decays into two excited states (^°'^Mop| and ^^'^Nd 6]). 

However, the calculation of the 2;//3/3-decay matrix elements is a complex task as the nuclei of interest are medium 
and heavy open shell ones. An additional complication arises from the fact that this is a second order process in 
weak interactions and, therefore, the construction of the complete set of states of intermediate nucleus is needed. 
At present there are basically two different approaches to the calculation of the /35-decay matrix elements which are 
the Quasiparticle Random Phase Approximation (QRPA) and the shell model j^l^. Both approaches have their 
specific advantages and disadvantages. In particular, the shell model is restricted to configurations with low-lying 
excitations, but takes into account all the possible correlations. On the other hand, the model space of the QRPA 
includes orbits far away from the Fermi surface, what allows one to describe high-lying excited states (up to 20 or 30 
MeV) of the intermediate nucleus. However, it is questionable whether this approach takes into account the relevant 
ground state correlations in a proper way. Also, the 2i^/3/3-decay matrix element evaluated within the QRPA exhibits 
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strong sensitivity to the details of nuclear Hamiltonian and dependence upon the model assumptions. The shell model 
has been found successful in the case of the 2^/3/?-decay of ^®Ca 7]. However, for other nuclear systems the shell 
model predictions are not convincing enough. In particular, for heavier nuclei it is not clear whether the considered 
model space is sufficiently large. This question becomes irrelevant only in the case when the dominant contribution to 
the 2vPf3-decay matrix element comes from the transitions through the lowest states of intermediate nucleus. Good 
candidates with this property are especially those nuclear systems for which the ground state of the intermediate 
nucleus is 1+ state, e.g., ^^°Mo, "^Cd and '^'^^Te. 

Some times ago it was suggested by Abad et al. that the 2i//3/3-decays with 1+ ground state of the intermediate 
nucleus are solely determined by the two virtual beta decay transitions: i) the first one connecting the ground 
state of the initial nucleus with 1^ intermediate state; ii) the second one proceeding from the 1^ state to the final 
ground state. This assumption is known as the single state dominance (SSD) hypothesis. The main advantage of the 
phenomenological approach based on the SSD hypothesis is that it is model independent since the values of the single 
/3-decay matrix elements can be deduced from the measured log ft values. 

The available experimental data0, 0, Ull IT^ and the results of theoretical analysis [T^l IT^ support the 
conjecture that the SSD hypothesis is valid for two-neutrino accompanied double beta decay. Recently, it has been 
suggested that the SSD hypothesis can be confirmed or ruled out by the precision measurement of differential rates 
in the case of 2j//3/3-decay of ^°°Mo The corresponding effects have been studied for the NEMO 3 experiment. 

Previous theoretical studies 'l4l of the SSD predictions were based on the approximation, when the sum of two 
lepton energies in the denominator of the 2i^/?/3-decay matrix element is replaced with their average value. In Rcf. 
UM it was shown that this treatment can be inadequate for some nuclear systems and demonstrated how to calculate 
the matrix elements without this approximation. 

In this article we present the results of systematic study of the neutrino emitting modes of double beta decay within 
the SSD approach without any additional approximation. The calculated decay rates are given both for transitions 
to the ground and excited final states. We also examine the question which of double beta decaying nuclei have the 
best prospects for experimental verification of the SSD hypothesis via the analysis of the differential rates. 

In Section |n] we describe our SSD based approach. The double beta decay rates are derived in Section ITTll Results 
and discussions are presented in Section ITVI Conclusions are drawn in Section Ivl 



II. SSD APPROACH TO DOUBLE BETA DECAY 



We distinguish four types of neutrino (antineutrino) accompanied nuclear double beta decays (/3/?-decays). The 
most attention has been previously paid to the mode with emission of two electrons (2j//3~/3~-decay), 

{Z~2,A)^ {Z,A) + 2e- +2i>e, (1) 

where Z and A denote atomic and mass numbers of a daughter nucleus, respectively. Another possible mode is double 
beta decay with emission of two positrons (2i//3+/3+-decay), 

{Z + 2,A)^{Z,A) + 2e+ + 2i^e. (2) 

In comparison with the 2zy/3~/3~ -decay this process is disfavored by smaller available kinetic energy and also by the 
Coulomb repulsion on the positron. More favorable for detection than the 2i'/3+/3+-decay are the capture processes, 
namely the capture of one bound atomic electron with emission of positron {2i'ef3~^ -decay) 

e- + {Z + 2,A)^ {Z,A) + e+ + 2v,, (3) 

and the double capture of two bound atomic electrons (2;/ee-decay) 

2e^ + {Z + 2,A)^ {Z,A) + 2v,. (4) 

The enhancement of the capture processes in comparison with the 2v (3^ (3^ -decay originates from the larger available 
energies and from the fact that the suppression due to the Coulomb repulsion is partially (2zye/3+-decay) or fully 
(2£e-decay) avoided. The electron capture from the s-states is a dominant subprocess, because the corresponding 
electron wavefunction does not vanish at the origin. In practice, one takes into account the capture both from the K 
and Li shells. However, the probability of capture from Lj shell is suppressed by about a factor 10 in comparison 
with capture from K shell. The inclusion of the electron capture processes into our analysis is, in particular, justified 
by a reviving interest to the experimental study of these processes ^3 • 
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The nuclear double beta decay is a second-order process of perturbation theory in weak interaction. The standard 
/9-decay Hamiltonian has the form 



^P{x) = -^e(:r)7^(l - l,)ve{x)j^{x) + h.c. 



(5) 



where Gp — Gp cos 0c and Gp is Fermi constant, 9c is Cabbibo angle, e[x), i'{x) are respectively electron and neutrino 
field operators, j^{x) is strangeness-conserving free charged hadron current. 

The /3/?-decay matrix elements involve summation over the discrete states and integration over the continuum states 
of the intermediate nucleus in the presence of energy denominators corresponding to the second order of perturbation 
theory in weak interactions. The neutrino accompanied modes of /3/3-decay are dominated by the Gamow- Teller 
transitions allowed by the angular-momentum, parity and isospin selection rules. The contribution of double Fermi 
transition to the double beta decay amplitude has been found to be small As is known, in nuclear models, like 
shell model and QRPA, the evaluation of the nuclear matrix elements is restricted to the two subsequent Gamow- Teller 
transitions through a certain number of low-lying discrete states of the intermediate nucleus. 

The partial /3/?-decay matrix elements associated with the transition to the ground or excited states of the final 
nucleus with angular momentum and parity J'^ =0^, 2+ can be written as 



3 ^F'^^ 



IIE,'TSa,. II 0+) 



U,k = 1,2) 



(6) 



with s = 1 for = 0+ and s = 3 for J'^ = 2+. |0^), \JJ) and |1+) are respectively the wave functions of parent, 
daughter and intermediate nuclei with the corresponding energies Si, ej and £„i- The terms Ej (j = 1,2) represent 
the energy of outgoing electron/positron or (with the minus sign) the energy of bound electron and Vk {k = 1,2) is 
the energy of emitted neutrino or antineutrino. Both, nuclear and lepton energies are assumed to be in the units of 
electron mass rrie- 

For the sake of simplicity the sum of outgoing lepton energies in the denominators of the matrix elements is 
usually replaced with its average value fwo, where / = 1/2 for 2v(3^(3^ and 2!^£e-modes and / = 1/3 oi f — 2/3 
for 2i/e/3+-mode {wq = {si - £/) = Wo/nie, see TablePl. This approximation allows one to factorize the lepton and 
nuclear parts in the calculation of the /?/3-decay rate. It works generally quite good for the cases when the dominant 
contribution to the /3/3 matrix elements comes from the transitions through the higher lying states of the intermediate 
nucleus. 

The SSD hypothesis drastically simplifies the calculation of the two-neutrino f3f3 nuclear matrix elements. It is 
supposed that nuclear matrix element is governed by the two virtual transitions: the first one going from the initial 
ground state to the 1+ ground state of the intermediate nucleus and second one from this 1^ state to the final = 0+ 
or 2+ state. With this assumption the relevant nuclear matrix element is as follows: 



M, 



(±) 

GT,m=g.s 



E 



(i: 



E 

j' 



(7) 



The superscript — (+) indicates the neutrino (antineutrino) emitting mode of the P(3 decay. The value of this matrix 
element can be determined in a model independent way from the single /3 decay and electron capture measurements. 
^From the experimental values of log ft for electron capture and single /3 decay of the ground state of the intermediate 
nucleus with J" = 1+ we get 




\- - lln+\ 1 '-^^ 



j' 



1 3D 



9 A V ftp- ' 



(8a) 
(8b) 



where D = (37r'^ln2)/(G^mg) is beta decay constant. The /3/3-emitters, which fulfil the SSD hypothesis condition, 
namely those with J"^ = 1+ ground state of the intermediate nucleus, are listed in Tables III and HTl 

In the previous SSD studies of neutrino (antineutrino) /9/3-decay the approximate treatment of the energy denom- 
inators has been applied [l^ ll^ . The sum of lepton energies entering the energy denominators were replaced with 
a half of the energy release of this process. Recently, it was shown that this approximation may lead to a significant 
inaccuracy in determining /3/3-decay half-lives for the transition to ground and excited states [T^ . It was demon- 
strated that the exact treatment of the energy denominators is especially important for /3/3-nuclear systems with a 
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TABLE I: Basic characteristics of P~ (3~ unstable nuclei, with 1"^ spin-parity of the ground state of intermediate nucleus. 
Wo — Ei — Ef is the energy difference between initial and final ground states in MeV deduced from atomic masses |2^ . 
A = El — Ei with El being the energy of the ground state of intermediate nucleus. 5 is the natural abundance of the isotope. 



™Zn «°Se iooMo 


^°4Ru 


llOpj 


"^Cd 


"«Cd 


Te 


Wo [MeV] 2.0229 1.1559 4.0563 
A [MeV] 0.1437 1.3596 -0.3429 
S 0.62 49.61 9.63 


2.3216 
0.6302 
18.62 


3.0217 
0.3814 
11.72 


1.5588 
0.9409 
28.73 


3.8270 

-0.0410 

17.49 


1.8892 
0.7408 
31.74 



TABLE IL Basic characteristics of experimentally interesting double electron capture unstable nuclei with 1^ spin-parity of 
the ground state of intermediate nucleus. E^{K) and Ee{L) are the energies of bound electron for K (lsi/2) and Lj (2si/2) 
shells, respectively. ?lo,-i (%,-i) denotes the probability to find lsi/2 (2S1/2) state electron in the volume around the 

origin of nucleus (see Eqs. 11231 and II311 '1. Wo, A, Ee{K) and Ee{L) are given in units of MeV. The quantities Wo, A, 5 are 
defined in Table 13 





"Zn 


™Kr 


"'«Cd 


io«Cd 


"^Sn 


i20Te 


i^OBa 


Ce 




">*Er 


Wo 


0.0743 


1.8440 


1.7491 


-0.7528 


0.9002 


0.6764 


1.5886 


1.3752 


0.8225 


-0.9979 


A 


1.0897 


1.2188 


0.7051 


2.1600 


1.1747 


1.4932 


0.8796 


0.9838 


0.8067 


1.4738 


5 


48.63 


0.35 


1.25 


0.89 


0.97 


0.09 


0.106 


0.185 


0.14 


1.61 


E4K) 


0.4986 


0.4931 


0.4786 


0.4786 


0.4758 


0.4728 


0.4664 


0.4630 


0.4437 


0.4437 


E^L) 


0.5079 


0.5065 


0.5028 


0.5028 


0.5021 


0.5014 


0.4997 


0.4988 


0.4939 


0.4939 




4.521 10~^ 


8.792 10"^ 


2.776 10"2 


2.774 10"^ 


3.311 10~2 


3.934 10"^ 


5.524 10~2 


6.522 10"^ 


1.467 10"! 


1.465 10-1 


^1,-1 


5.896 10"^ 


1.168 10"^ 


3.872 10"^ 


3.869 10"^ 


4.662 10"^ 


5.592 10"^ 


8.020 10"^ 


9.574 10-3 


2.291 10-2 


2.288 10-2 



small energy difference between the intermediate and initial ground states in comparison with the value of Wq, e.g., 
A=100 and 116 systems. 

The effect of energy denominators may have a strong impact on the differential rates. In Ref. |1(t| [ the two limiting 
cases have been considered: i) The differential rates have been evaluated within the SSD hypothesis with the exact 
treatment of the energy denominators, which are the functions of lepton energies, ii) The above mentioned standard 
approximation for the energy denominators has been applied in which case the energy denominators are independent 
of lepton energies. 

As we commented above, this approximation is justified if the transitions through the higher lying states of the 
intermediate nucleus give the dominant contribution to the /3/3-decay amplitude. This assumption, which we call as the 
higher state dominance (HSD) hypothesis, is an alternative to the SSD hypothesis leading to the different predictions 
for the differential rates. Thus, from the comparison of experimental data on differential rates with SSD and HSD 
predictions one can confirm or rule out the SSD hypothesis. This would provide us with a valuable information about 
the structure of involved nuclear matrix elements. 



III. DOUBLE BETA DECAY RATES 

In what follows we present basic formulas needed for evaluation of the /3/9-decay half-lives and the differential rates 
within the SSD and HSD hypotheses. 

A. The 2!^/3*/3*-decay half-life 

The differential rate for 2vl3~ and 2z//3+/3+-decays leading to (J^ — 0+, 2+) state of the daughter nucleus 
takes the form [Ull^j 

dujff^'^^ = a2uA%'^^'^^ {ei,e2,vi,V2)dVt, (9) 
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where 021/ = (G'^(7A)^'7ig/(647r^). The phase-space factor dfJ is 

dil = 2nin2eie2i^fi'iS{ei + 62 + + V2 ~ Ei + ef)deide2dvidv2. (10) 

Here, tt^ is the momentum of k-th emitted beta particle. The electron/positron momentmn ttj., energy and neutrino 
energy (fc = 1, 2) are assumed to be in the units of electron mass m^. The function y^j^ can be written as 

A^jKP {61,62,1^1,1^2) ^ a{ei,e2)Mff ^ (61,62,1^1,1^2) (11) 

with 

a{6i,62) = )(£i)4f (£i)^^o^*^(£2)i?i?(e2), (12) 

Here, the superscript — (+) denotes 2i^/3~/3~ (2i^/3+/3+) - decay mode. i^Q ''(e) is the relativistic Coulomb factor [Tgj . 
The correction factor r[-^^\6) is given by 

r[-\6) ^ R-l{6) + R+i{6) (13) 

with i?±j,(e) = Z?j_j.(e)(e =p TOe)/(2e). The function Z?^j,(e) reflects the charge distribution inside the nucleus. Its 

explicit form can be found in Ref. The calculation of Fq'^\6) and R^\6) differs from the calculation of Fq •*(£) 

and r[i ''(e) in the sign of Z associated with the final nucleus. 

The function Ai^^J^ ^ {61, 62, vi, 1^2) in Eq- lll|) is the product of the lepton energy dependent nuclear matrix 
elements. Within the SSD hypothesis it takes the form 

M'fJ'^''^{6i,62,Vl,V2) = |M^^)g.,.(J'^)|2/Cj.(£i,e2,I^l,I^2). (14) 

The nuclear matrix element ^ ^ [J'^) consists of two single /? decay matrix elements associated with transitions 
via 1+ ground state of the intermediate nucleus [see Eq. JTJ]. The energy denominators enter into the expression for 



A:j-(£i,e2,t^i,^2) = -{K^,s.{^i, 62, 1^1, 1^2) + LI s,{6i, 62,1^1, 1^2) 



1 

+Kg,sX£i, £2, fi,i^2)Lg.s.{6i, 62, 1^1,1^2)) {for J'' = 0+), 
(i^g.s.(£i,£2,i'i,i^2) - Fg.s.(ei,£2,J^i,i^2))^ {for J"" = 2+) (15) 



with 



Km=g.s.{£l^£2,l^l,1^2) = {Sm^g.s. - £i + £l + l^l) ^ + {Sm^g.s. - £i + £2 + ^2) \ (16) 
Lrn=g.s.{£l,£2,l^l,1^2) = {£ni=g.s. - £i + £l + 1^2)^^ + {£rn=g.s. - £i + £2 + Vl)^^ ■ (17) 

Here, 6m=g.s. corresponds to the energy of the 1"*" ground state of the intermediate nucleus. 

Let us note that the general matrix element used in the nuclear model dependent calculation of the 2z^/3^/3^- 
decay rates includes the sum over the intermediate nuclear states. In this case it is evaluated within the standard 
approximation, when in the common denominator of Km and Lm the lepton energies are replaced with their average 
values 

(2i//3±/3±) ^ ^ 

61 ~ 62 — Vl — V2 — , Wg — £i — £/, (_l5j 

for the transitions to the 0+ final states. For the transitions to the 2+ final states the above replacement is applied 
to the lepton energies in the denominator of the expression 

^ __j _ 2(t/i - V2){6i - £2){£m -£i + {£i - £/)/2) , . 

(£m - £i + £1 + l^l){£m ~ £i+£2+ V2){£m - £i + £l + i^2)(£m - £i + £2 + Vl) 

As the result of these approximations the expression for M^j'f ^ in Eq. 114|l can be factorized into the product of 
the lepton-energy independent sum over the intermediate nuclear states and the factor that does not depends on £,„. 
Its explicit form for J'^ = 0+, 2+ is given in Table UTTI 
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TABLE III: General nuclear matrix elements Mjtt entering into the expressions for a — 2i//3^/3*, 2uef3^ , 2!/££-decay rates [see 
Eqs. H25|l and l|8()|l ] with the standard approximation in their common denominators. For the explanation and definition 

of Wg"' see Eqs. 1181 . (1281 and 1341 . The quantities ~ £m — £i include the energy of intermediate l""" nuclear state 
and the ground state energy Si of initial nucleus. The summation "^^m runs over the intermediate nuclear states. Sei with 
I = 1,2 stands for the energy of captured electron, ei and Vk with l,k = 1,2 are the energies of emitted electron/positron and 
(anti)neutrino, respectively. All the energies are in the units of rrie. 



mode J" = 0^ 



J" = 2^ 



2i^/3±/3^ 



2i/e/3^ 



2i/£e 



A 

E™mW„(o+) 



+ 



+ 



— 

^2 2 



T MyT> (2 



{±) 1^0+ \ 2(ei -E2)('^l-'^2) 



A„ + - 



{+) ^^,+ ^ 2(ei-£2)(t'l-'^2)(A^ + (£,-ej)/2) 



A™ + ^ 



,,{+) 2{ee2-eei)(i'l-"^2){A„ + {ei-ej:)/2) 



„(<») ^ 



A„~e„+- 



Within the SSD hypothesis the total 2i//?*/?*-decay rate can be written as 

uiff ^ ■)„f!±r± , — : = cf^'" / deiEiTTi / d£2£27r2 / dviv(v2a{ei,e2)lCj^{£i,e2,vi,V2), (20) 



-'SSD W 

± -L U 

where a = 2i//3^/3='=. The factor cf^^ determined from Eqs. Q and ||HJ) takes the form 

^SSD _ 9„ \^j± |2 _ 9 (?n2)^ 

C2^ - ^a2i,\M(jj, \ - — ■ (21) 

It is worthwhile to notice that the f3^ -decay half-life in Eqs. (|2()|l and H21|l does not depend explicitly on the 
values of and the axial weak coupling constant gA- 



B. The 2!/e/3+-decay half-life 

In the similar fashion we treat 2;ye/3+-decay mode with 0+ and 2+ states of final nucleus. In this case the differential 
decay rate can be written in the form |2Clj| 

diufj^^ = a2,N^'''^^ A^'^^-Se,, 81,1^1, i^2)dn. (22) 

The factor N'^^'^^* is associated with the normalization of wavefunction of bound electron, namely 

^2ue!}+ ^ f 2(27r2)01o -1 K{lsii2) capture 
\ 2(27r2)01i _i L/(2si/2) capture 

Here, 01o,-i (^i,-i) is the probability to find lsi/2 (2si/2) state electron in the volume l/m^ around the origin of 
nucleus. The details of the calculation of these factors can be found in |23|. In Table HTl we present their numerical 
values for some nuclei. One may note that for L/-shell the above probability is suppressed by about a factor of 1/8 
in comparison with that for iiT-shell |23l- The 2r/e/3+-phase-space factor is given by 



dfi = 2'Kieiu1v25{ei — + vi + 1^2 ~ £i + e f )deidvidv2, 



(24) 
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where £ei is the energy brought to nucleus by a captured electron. Next, we have 

A%'''^^ {-Ee:,, s 1, 1^1,1^2) ^ a{ei)M%'"^^ {-Ee,, 81,1^1,1^2) (25) 

with 

a(£i)=i^(+)(ei)i?W(£i). (26) 

The SSD form of the matrix element is 

M%'^\-ee„ 61,1^1,1^2) = \M^GT,^ jr)\^ICj.{ei,-ee„>^i,>^2), (27) 

where the nuclear matrix element and the energy denominator in the r.h.s. of this equation are defined in Eqs. ((7|) 
and H15|l. respectively. 

For a comparison with the SSD expression for M'^^i^^ in Eq. (|27(l we present its standard form in Table ITTll which 
includes transitions through the excited states of intermediate nucleus and is based on the conventional approximation 
for the lepton energies in the energy denominators 

(2i/e/3+) 

£1 ~ I/l ~ Z/2 — ^^-^ , '=£i-£/+eei. (28) 

Returning to Eq. (|22() . we write down the total 2j/e/3+-decay rate within the SSD hypothesis as 



„(-) 



2i.e/3+ lri{2) 

(JJ IT, 



j^2uep+ (•jit -J 

where a = 2vej3^ . 



— =cl^''N^'''f^^ J deiSiTTi J dvivlvla{ei)lCj.{ei,-e,„v^,v2), (29) 



C. The 2i/ee-decay half-life 

The double electron capture differential rate is given by the formula poj 

dw2.se ^ a2uN^''''M%''{-ee„-ee„vi,v2)dn (30) 

where 



T2uee _ j (27r2)2m2 _i KK Capture 

2(27r2)2fno KLj capture 



j^2uee ^ ^ V^" ^ -"^0,-1 " Ldpbuie ^^^^^^ 



= 2v1v25{—ee^ — £e2 + ^1 + ^2 — £j + e f)dvidv2- (32) 

The probability factors 01o,-i, have been defined in the previous subsection. Their numerical values for some 

nuclei are shown in Table ITU 

Within the SSD hypothesis the matrix element j';^'^^(— , —Se-ii vi, 1^2) in Eq- (|30|l takes the form 

Mf/%-ee„-ee„Pi,iy2) = M^^^^^^{J^)lCj.{-ee„-ee„vi,v2). (33) 

Here, £ei , ^^'^ the energies of bound electrons captured from the K or i/ atomic shell. The nuclear matrix clement 
and the energy denominator in the r.h.s. of this equation are defined in Eqs. (U)) and IjlSI) . respectively. 

In the standard treatment of the 2i^£e-decay the nuclear matrix element (see in Table ^QJ includes the contribution 
from transitions through the excited states of intermediate nucleus, which are constructed within an appropriate 
nuclear model. It is also assumed that 

{2vee) 

Vi~V2~^^ , ' = £i - £/ + £ei + £e2. (34) 
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This approximation allows one to separate the phase space integration from the evaluation of the 2i/e£-decay nuclear 
matrix clement. 

In our derivation of the SSD predictions we always take into account the dependence of the corresponding nuclear 
matrix elements on the lepton energies as in Eq. Thus, the SSD expression for the total 2:/ee-decay rate derived 

from Eqs. ^U^, reads 



^ssD^2.se J dv^vlvlKj.{-ee,,-ee,,i^i,v2) (35) 



with a = 2vse. 



ln{2) 



D. Energy distributions of outgoing electrons 

There is a unique possibility to prove or rule out the SSD hypothesis for /3/3-decays by experimental measurement 
of the energy distributions of outgoing electrons 16]. Any deviation from the SSD prediction for the shape of the 
electron energy spectra can be interpreted as the influence of the transitions through excited states of intermediate 
nucleus. The HSD hypothesis, mentioned in Section|nl reflects an opposite situation to the SSD one. It corresponds 
to the approximation when one neglects the dependence of the relevant nuclear matrix element on the lepton energies 
which is well justified if the dominant contributions to the /?/3-decay rate come from virtual single-/3-decay transitions 
via the higher lying 1+ states of intermediate nucleus. In this case one substitutes the lepton energies in the energy 
denominators of nuclear matrix elements with their average values defined in a special way as in Eqs. ((TH|l . 

For the experimental verification of the SSD and HSD hypotheses the most appropriate nuclear systems are those, 
which exhibit a maximal difference in the shapes of the electron/positron energy spectra predicted by the SSD and 
HSD based scenarios, meanwhile having relatively low /3/3-decay half-lives. 

Thus, the subject of our interest are the differential rates normalized to the total decay rates, namely 

Pj"^''(£:) = (/ = SSD, HSD). (36) 

For 2i^/3^/9^-decay {a = 2vf5^f3^) £ stands for the single electron/positron energy £i or for the sum of the kinetic 
energies of outgoing electrons/positrons t. For 2:/e/3+-decay (a = 2ve(3^) £ denotes the positron energy ex. It 

is worthwhile to notice that the quantity Pj^^'^ {£) does not include any nuclear matrix element because both the 
differential and total rates are scaled with the same nuclear matrix element which cancels out in their ratio. The only 
nuclear structure inputs are the energy release in the considered process and the energy difference between 1"*" ground 
state of intermediate nucleus and parent nucleus. 

We proceed with presenting the expressions for the single electron distribution dujjT,/dei and the summed electron 
spectrum duij^/dr for different /3/3-decay modes within the SSD and HSD hypotheses. We consider experimentally 
most favored ground-state transitions. 

The 2v(3~ j3~ -decay to ground state: 
i) The single electron differential rate 



dsi 



2a2,.|Mgj, g^ pEiTTi J de2e2TT2 J dvivlvla{ei,e2)K,Q+ ^ {ei, 62, 1^1,^2) (SSD), (37) 
1 

u)o°'-ei u)o°'-ei-e2 

= 2a2,|7W^|°pei7ri J de2e2^2 j di^ii^fi^la{ei, 82) (HSD). (38) 
1 

Here a = 2vj3~ I3~ . We stress again that the corresponding SSD and HSD total decay rates are proportional to the 
squared nuclear matrix elements Mqj,^^ and M^^^ , respectively. Therefore, they drop out from the normalized 

differential rates Pj"^'^(ei) in Eq. 
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ii) The summed electron differential rate 

" ^ 2a2^|MQj,g^ P / deiEiTTi / dz^iz^f i/|a(ei, e2)/Co+^ (ei, £2, i^i, 1^2), {SSD), (39) 



dT 



1 

r+l u)^,°'-r-2 



= 2a2^|M^|f|^ y deieiTTi y dvxvlvla{ei,e2) {HSD). (40) 
1 

Here, a = 2vj3^ fi^ and r = ei + £2 ~ 2 is the sum of the kinetic energy of emitted electrons in unit of mg. 

The 2vel3^ -decay: 
The single positron differential rate 



.„(=) 



-El 

duj 



= 2a2.|M+j,g^ piVei^i / di/izy2^|a(ei)/Cj. (ei, -Ee^ , z^i, 1^2) (5^^?), (41) 



= 2a2^iV|7W|fFpei7ri J dvivlulaisi) {HSD). (42) 



Here a = 2z^e/3+. 

In this paper we do not discuss the single positron energy distribution and summed positron energy spectrum of 
the 2i//3+/3+-decay. As we commented in Sectional this process is less favorable for experimental study than the 
other /3/9-decay modes because of its small available energy and the Coulomb suppression. Also we do not discuss the 
angular distributions in 2i^/3^/3^-decay which may also be helpful for the experimental verification of the SSD 
and HSD hypotheses. 

There may happen the situation when the measured energy distributions do not correspond neither to the SSD nor 
HSD predictions, so that the SSD is only partially realized and there are interfering contributions from the higher- 
lying 1+ states of intermediate nucleus. In this case one may apply to the data analysis the approximate expression 
for the mixed SSD-HSD distribution: 

pj:^'""^-(f)«.TPje^'^^^(f) + (l-a;) Fje^'^^''(f) (0<a;<l), (43) 

where the parameter x characterizes the relative contribution from the first 1+ state of the double-odd nucleus to 
the corresponding mode of /3/9-decay. Fitting the above formula to the experimental distributions one can extract 
the value of x and get an idea on the role of this contributions in the studied process. However, one should take into 
account that the accurate definition of the mixed distributions in Eq. H43|l requires the knowledge of nuclear matrix 
elements involved in the differential rates. 



IV. RESULTS AND DISCUSSION 



In the present work we have analyzed the neutrino accompanied double beta decay modes for emitters, which fulfil 
the SSD hypothesis condition: The ground state of odd-odd nucleus must be 1+ state. The subject of our interest are 
eight even-even nuclei, which satisfy this condition for the 2zy/3~/3~ -decay, and ten double electron capture unstable 
nuclei. Their basic parameters are listed in Tables HI and Hll iFrom Table HI one may notice that in the case of ^'^'^Mo 
and ^^^Cd the energy difference A between the ground states of initial and intermediate nuclei is negative. Therefore 
for these nuclear systems the difference between the maximal and minimal values of the energy denominators of 
/3/3- nuclear matrix elements is significant. In this case it is especially important to properly take into account the 
dependence of nuclear matrix elements on lepton energies in the calculation of decay rate. 

In Table HTl we listed the dimensionless normalization constants 01o,-i a-nd ^i.-i 20] of relativistic wave function 
of bound electron in the K (lsi/2) and Lj (2si/2) states, respectively. Comparing the values of these constants one 
can see that the electron capture from Lj shell is suppressed by about a factor of 10 in comparison with the electron 
capture from K shell. The energies of bound electron and normalization quantities 01o,-i and Oli _i were calculated 
assuming the Coulomb potential for finite-size spherical nucleus (R — vqA^^^, tq = 1.2 fm) with uniform charge 
distribution. 



10 



TABLE IV; The effect of the correction factor R^j(e) (see Eq. [HJ on the calculation of 2///3*/3*-decay rate. The difference 
between the results with and without the correction factor is given in percents. Only transitions to the ground state of final 
nucleus are considered. 



™Zn «°Se 


i°»Mo 


i°*Ru 


llOpj 


ii^Cd i«Cd 




^»Kr i»^Cd 


i3°Ba "6ce 


2.6 3.5 


4.6 


5.7 


6.0 


7.0 6.2 


8.1 


19 33 


42 42 



We calculated the half- lives for the two- neutrino double beta decays of the nuclei listed in Tables HI and ^ within 
the SSD hypothesis. For this purpose we used the /i-values corresponding to the single f3~ and electron-capture 
transitions feeding the low-lying states of initial and final nuclei. They were taken from NNDC On Line Data Service 
of ENSDF database [l^. For those nuclear systems, where only the lower limits on the logft-value is known, we 
derived the lower limits on the /9/3-decay half-lives. The calculation of /3/3-decay rates has been performed with the 
exact treatment of the energy denominators of perturbation theory. Transitions to the final 0+^ ground and 0+ and 
2+ excited states have been considered. In the case of 2i/£:/3+-decay and 2r/ee-decay modes the captures from K and 
Lj shells have been taken into account. We neglected the captures from Ljj (2^1/2) and Lju (27)3/2) shells suppressed 
by the factor {r/R) in comparison with s-states. 

We improved our previous results [islIT^ by more accurate treatment of the Coulomb correction factor R^i{s) [see 
Eqs. H12|l and (|13|l ]. which reflects the charge distribution of the final nucleus. We have found that the presence of this 
factor affects the 2i//3~/3~ -decay half-life only slightly by reducing its value up to 5— 10%. On the other hand, it plays 
an important role in the evaluation of the 2i'/3+/3+-decay half-lives enhancing their values by 20 — 50% (Tab. llV)l . 

The results for the 2i//3"/3~-transitions to the final 0^^ ground and 0^, 2+ excited states are shown in Table Ivl 
It is seen that in the case of ^^^Mo the SSD hypothesis predictions are consistent with the experimental half-lives 
both for the transitions to ground and 0^ excited states. On the other hand for ^^^Cd the calculated half-life to the 
final ground state is by about a factor 2 shorter than the experimental value. A similar situation occurs also for the 
2v I3~ j3~ -dececy of ^^^Te. One possible explanation of this discrepancy may consist in the cancellation between the 
transitions passing through the 1"*" ground state of double odd nuclei and the rest of the transitions going through 
the excited 1"*" states. However, one should also keep in mind the uncertainties both in the measured 2z//3~/3~ -decay 
half-life and in the experimental determination of the log ft value of electron capture. In view of these uncertainties 
it is not possible to make definite conclusion on the validity of the SSD hypothesis for the case of this nucleus. For 
other 2z//3~/3~ -transitions in Table IVl there are only lower limits on the half-life available. 

In these cases the SSD hypothesis helps us to give an order of magnitude estimates for the the 2i//3~/3~ -decay 
half- life. One may note that the 2z//3~/3~ -decays of -^^^Mo and ^^^Cd to excited state are significantly faster than 
those to the 2^ state in spite of the smaller energy release. The suppression of the 2+ decay channel is due to the 
mutual cancellation of K and L terms in the 2v I3~ j3~ -decay amplitude [see Eq. (|15|l ]. 

The results for the 2i//3+/3+-, 2ve(3^- and 2i'££-decays are listed in Tables I VII IVIII and IVIIII So far none of these 
/3/3-decay modes have been seen experimentally. The SSD hypothesis calculations of decay rates require the knowledge 
of the measured Zoiy/i-values for single /3~ and electron-capture decays. However, the /of;/i-values for the /3~ decay of 
the ground state of odd-odd nucleus are not known for A — 120, 136, 162 systems and, therefore, we do not present the 
SSD predictions for the half-lives of the corresponding processes. In the case of ^^''Er its small available energy allows 
only 2ve£ mode, which is strongly suppressed. Therefore, the SSD predictions for this transition are not presented. 
In particular cases, only lower limits on the log ft value for the single /3~-decay of the ground state of intermediate 
nucleus are measured. This allows us to determine only the lower limits on the corresponding half-lives within the 
SSD based approach. At present the attention of experimentalists is concentrated on the 2i^ee-decay of ^°^C(i jl^ . 
For the ground state to ground state transition of this process the SSD hypothesis predicts the half-life larger than 
4.4 X 10^^ years. 

In Fig. n we present the single electron differential decay rate Pj"!)'^{£) [see Eqs. (|36|l . (|37|l and (|37|l ] for the 
2v I3~ j3~ -decay of ^"^Afo, ^^'^Pd and ^^^AIo to the 0+ ground state. We see that there is a notable difference in the 
behavior of the single electron energy distribution calculated within the SSD and HSD approaches especially for the 
small electron energies. The effect is significant in the case of ^"°A/o and ^^^Cd. This is due to the smallness of 
the energy difference between the ground states of the intermediate and initial nuclei in comparison with the energy 
release in the 2i//3~/3~ -decay process (see Table IJ). We also note that the 2iy f3~ f3~ -decay of ^"^Afo is slightly more 
favored for the experimental verification of the SSD hypothesis than the 2z//3~/3~ -decay of ^^^Cd due about 3 times 
shorter half-life (see Table IVjl . 

The difference between the SSD and HSD predictions at the beginning of the single electron spectra is of the 
order of one percent. Thus the experimental observation of this effect requires large experimental statistics. The 
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TABLE V: The half-lives (rggnf ^ for the 2u f3~ /3~ -decay transitions to the ground state (J'^ — O^.s.) and excited states 
[J'" — 0'^ and 2"'") calculated within the SSD hypothesis. The excitation energies of the states, E'^^'^, and their log fti (electron 
capture) and log/t/ (single /3~-decay) are shown. The corresponding values were taken from NNDC On Line Data Service 

from ENSDF database updated on 08/13/2003 Q. Here Wo = m^w^^'"'^' '^'^ (see Eq. lO). The experimental 2z//3/?-decay 

half- lives rixp'' ' for ^°°Mo to ground and excited 0^ states are from Refs. |24| and f2^ . respectively. The 2z//3/3-decay 

half- life of ^^^Cd is from Ref. [2R| |. For the lower limit on rixp'' ' we used the best limit presented in Tables of Ref. 3, 
which is a compilation of experimental data. 
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TABLE VL The half-lives {T^^ ^ ') for the 2i//9^/3^-decay transitions to the ground state {J^ = 0^,,.) and excited states 
[J^ — 0"*" and 2"*") calculated within the SSD hypothesis, logj^. and logj>t^ denote iog/t- values for single and electron 

capture-decays, respectively. Here Wo ~ uieW^^'^ ^ ^ (see Eq. CHJ)- For other notations see TabfVl 
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0.0 


4.92 
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> 2.4 10^" 
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5.36 


2+ 
^1 


536.95 


6.3 


1.0516 1.2 10®° 


> 4.0 10^^ 



"Only the lower limits of half-lives are calculated 

corresponding conditions have been reached in the NEMO 3 experiment with ^^^Mo collected more than one hundred 
events of 2i^/3/3-decay of ^'^"Mo 24]. The obtained results are in favor of the SSD hypothesis for transition to the 
ground state. There is an additional possibility for a cross-check of this conclusion by analyzing the differential decay 
rate P{T) [see Eqs. (|36ll . (|39ll and H40|) ]. where T is the sum of the kinetic energies of outgoing electrons. For the 
2vl3(3-decay of ^^^Afo and ^^^Cd the summed energy distributions are shown in Fig. |21 We see that the shape of HSD 
and SSD curves is again quite different. The maximum of the HSD distribution is slightly above the maximum of the 
SSD one. The large statistics to be collected by the NEMO 3 experiment in the near future will allow us to study 
this effect as well. 

In Fig. O we show the single positron differential decay rate P{e) normalized to the total decay rate [see Eqs. 
ijSni and (gni] for ground state to ground state 2z^e/3+-decay transitions ^°'^Cd — > ^'"'Pd, ^^^Ba —> i^°Xe, and 
^^^Ce — > ^^^Ba. One can notice a significantly larger difference between SSD and HSD distributions in comparison 
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TABLE VII: The half-lives {Tgg'^^ ') for the 2!/e/3^-decay transitions to the ground state {J^ = O^.s.) and excited states 

[J^ = 0''' and 2+) calculated within the SSD hypothesis. Here Wq — m^Wi-^'^^^ ^ (see Eq. I28II '). The values in the first and 
the second rows for each transition correspond to the case of electron capture from the K and Li atomic shells, respectively. 
Other notations are the same as in TablVI 
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"Only the lower limits on half-lives are calculated 

with the case of the 2z//?~/?~-decay of ^"^Af o (see Fig. 0. This is partially due to the fact that in these two cases with 
emitted positron and emitted electrons the behavior of the relativistic Coulomb factor as a function of energy is quite 
different. In the case of positron its value is vanishing for kinetic energy equal to zero while for the electron it takes 
a maximal value. On the other hand the predicted SSD half-lives of the studied 2i/e/3+-decays (see Table IVIIII are 
about 2-3 orders of magnitude larger than the measured half-life of 2i//3~/3~ -decay of ^"^Mo. Thus, the experimental 
study of the 2:/£/3+-decays with the purpose of verification of the SSD hypothesis looks unfavorable. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have systematically studied the SSD hypothesis for all /3/3-unstable isotopes with the 1+ spin-parity 
ground state of intermediate nucleus presented in Tables HI Ull The two-neutrino accompanied /5~/3^, /3"^/3^, ECfi^ 
and EC /EC modes of double beta decay to ground (0+) and excited (O"*", 2+) states of final nucleus were considered. 
The half-lives of these processes were derived within the SSD hypothesis using the values of the matrix elements of 
single-/?^ and electron capture decays extracted from the measured log ft values. In those cases when only the lower 
limits on the log ft value were available, we calculated the lower limits on the corresponding double beta decay half- 
lives. The SSD calculations of /3/3-decay rates were performed without additional approximations keeping explic itly 
the dependence of energy denominators on the lepton energies |0| . In comparison with the previous studies 0, 0| 
we have more accurately taken into account the effect of Coulomb suppression in the part which corresponds to the 
charge distribution of final nucleus and described by the correction factor Rti{e) (see Eqs. 1121 126|) . We have found 
that the improved treatment of the Coulomb interaction of emitted charged lepton with the nucleus modifies the 
total rate of 2i//?~/?~-decay in 5 — 10%. For the case of the 2j//3+/3+-decay this improvement results in 20 — 50% 
enhancement of the half-life (Tab. IIV|) 

We demonstrated that the existing experimental data on the half-lives of the 2://3^/3~-decay of ^"^Afo (to O'^ 
ground and Of excited |5j|] states) and of ^^^Cd (to 0+ ground state 25]) show a clear tendency to favor the SSD 
hypothesis. However, we are not yet in the position to make definite conclusion on the validity of the SSD hypothesis. 
This is, partially, because of quite large uncertainties (about 50%) in the SSD predictions stemmed from insufficient 
precision of the existing experimental measurements of the logftEC values for electron capture. Thus, a more accurate 
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TABLE VIII: The half-lives {T^g'^^^) for the 2i'££-decay transitions to the ground state {J^ — Oj^.s.) and excited states (J'^ = 0"*" 
and 2"*") calculated within the SSD hypothesis. Here Wo ~ m^w'^^'^^^^ (see Eq. (I34ll 'l. For transitions to 0^ state the double 
capture of atomic electrons from K-shell (upper row) and from KLj (lower row) shells is considered. In the case of transitions 
to 2^ state only the double capture of atomic electrons from K-shell is taken into account. Other notations are the same as in 
Tab El 
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"Only the lower limits on half-lives are calculated 



experimental information on the associated /3- and /^/^-transitions is needed. 

The SSD approach allowed us to make the predictions for those neutrino (antineutrino) accompanied /3/3 decays to 
ground and excited final states, which have not yet been experimentally observed. Despite the fact that these predic- 
tions should be taken as order-of-magnitude estimates they may help in evaluation of the prospects of experimental 
searches for mentioned processes. 

The SSD results presented in Tables IVIVIIII show that there is a chance to observe 2z^ee-decay for ^"^Cd, ^^'^Sn 
and ^^°i?a at the level of 10^^ — 10^^ years. On the other hand for 2z^e/3+- and 2i//3+/3+-decays the SSD predicts, in 
general, significantly larger half-lives looking quite pessimistic from the view point of the experimental observation 
of these processes in the near future. We also found that with the exception of the 2:^/3" /3~ -decay of ^°°Mo to 0+ 
excited state the SSD half-lives of all the studied transitions to excited states 0+ and 1+ states are above 10^^ years. 

We proposed the study of the differential decay rates as a new perspective possibility for the experimental verification 
of the SSD hypothesis. One theoretically important point with respect to these characteristics consists in the fact that 
the differential rates normalized to the corresponding full decay rates do not depend on the values of the associated 
nuclear matrix elements. We have shown that the best candidates for this study are double beta decay chains with 
low- lying 1+ ground states of intermediate nucleus. In this case the shape of the single electron/positron distribution 
and the summed electron spectrum calculated within the SSD are very sensitive to the lepton energies present in the 
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FIG. 1: The single electron differential decay rate normalized to the total decay rate P{e) = {l/u)duj / de vs. the electron 
energy e for 2///3"/3"-decay to the 0+ ground state. The resuhs are presented for the cases of ^""Mo ^ ^""Ru, ""Pd "°Cd, 
and ""^Cd ^ "''Sn. The calculations have been performed within the single-state dominance hypothesis (SSD) and with the 
assumption of dominance of higher lying states (HSD). 



energy denominators. On the other hand, if the main contribution to the double beta decay matrix elements comes, 
contrary to the SSD hypothesis, from the transition through the higher lying states of intermediate nucleus (HSD 
hypothesis) the dependence on lepton energies in the energy denominators is negligible resulting in the distributions 
which can be experimentally distinguished from the SSD case. In the present paper we have shown that the precision 
measurements of the differential characteristics of the two-neutrino modes of the /3~/3~-decay of ^""Mo and ^^^Cd and 
EC/l3+-decay of i°^Cd, ^'^°Ba and ^'^^Ce are able to confirm or rule out the SSD hypothesis. Recently the NEMO 3 
collaboration have started the data analysis aimed to distinguish the SSD and HSD predictions for the 2v (3~ I3~ -Aec&y 
of i°°Mo [24]. 

We believe that this and further studies of the SSD hypothesis will allow us to improve the understanding of 
the structure of two-neutrino double beta decay matrix elements and will be helpful in fixing the nuclear structure 
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0.8 




T (MeV) 



FIG. 2: The differential decay rates normalized to the total decay rate P{T) = {l/uj)duj / dT vs. the sum of the kinetic energy 
of outgoing electrons T = me(ei + £2 — 2) for 2i//3/3-decay of ^°°Mo and ^^^Cd to the ground state of final nucleus. The 
conventions are the same as in Fig. 



parameter space needed for the calculations of the matrix elements of neutrinoless double beta decay within the most 
advanced nuclear models [271. 
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